import math
import numpy
import matplotlib.pyplot as plt
import scipy.fftpack
data = numpy.loadtxt("text.txt")
fs = 640e3
t = numpy.arange(0, (1 / fs * 12800), 1 / fs)
for index, value in enumerate(data):
        # if value > 127:
        #     data[index] = data[index] - 256
        data[index] += 2* math.sin( math.pi * 50e3 * t[index])

        data[index] += 2* math.sin(math.pi * 100e3 * t[index])
        data[index] += 2* math.sin(math.pi * 110e3 * t[index])

        data[index] += 2* math.sin(math.pi * 190e3 * t[index])
        data[index] += 2* math.sin(math.pi * 200e3 * t[index])
        data[index] += 2* math.sin(math.pi * 210e3 * t[index])

        data[index] += 2* math.sin(math.pi * 260e3 * t[index])
        data[index] += 2* math.sin(math.pi * 270e3 * t[index])
        data[index] += 2* math.sin(math.pi * 290e3 * t[index])
        data[index] += 2* math.sin(math.pi * 300e3 * t[index])

        data[index] += 2* math.sin(math.pi * 341e3 * t[index])
        data[index] += 2* math.sin(math.pi * 350e3 * t[index])
        data[index] += 2* math.sin(math.pi * 360e3 * t[index])
        data[index] += 2* math.sin(math.pi * 370e3 * t[index])
        data[index] += 2* math.sin(math.pi * 380e3 * t[index])

        data[index] += 2* math.sin(math.pi * 410e3 * t[index])
        data[index] += 2* math.sin(math.pi * 420e3 * t[index])
        data[index] += 2* math.sin(math.pi * 430e3 * t[index])
        data[index] += 2* math.sin(math.pi * 450e3 * t[index])
        data[index] += 2* math.sin(math.pi * 460e3 * t[index])
        data[index] += 2* math.sin(math.pi * 470e3 * t[index]) 

        data[index] += 2* math.sin(math.pi * 490e3 * t[index])
        data[index] += 2* math.sin(math.pi * 500e3 * t[index])
        data[index] += 2* math.sin(math.pi * 510e3 * t[index])
        data[index] += 2* math.sin(math.pi * 520e3 * t[index])
        data[index] += 2* math.sin(math.pi * 530e3 * t[index])
        data[index] += 2* math.sin(math.pi * 540e3 * t[index])
        data[index] += 2* math.sin(math.pi * 550e3 * t[index])

        data[index] += 2* math.sin(math.pi * 570e3 * t[index])
        data[index] += 2* math.sin(math.pi * 580e3 * t[index])
        data[index] += 2* math.sin(math.pi * 590e3 * t[index])
        data[index] += 2* math.sin(math.pi * 600e3 * t[index])
        data[index] += 2* math.sin(math.pi * 605e3 * t[index])
        data[index] += 2* math.sin(math.pi * 610e3 * t[index])
        data[index] += 2* math.sin(math.pi * 620e3 * t[index])
        data[index] += 2* math.sin(math.pi * 630e3 * t[index])

numpy.savetxt("data_rfi6.txt", data)
xn = scipy.fftpack.fft(data)
t = numpy.arange(0, (fs + fs/12799), (fs/12799)) + 1182e6
plt.plot(t, abs(xn/12800))
plt.show()